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Abstract. We present the analysis of the first two years 
of a search for microlensing of stars in the Small Mag- 
ellanic Cloud with the EROS (Experience de Recherche 
d'Objets Sombres) project. A single event is detected, al- 
ready present in the first year analysis. This low event 
rate allows us to put new constraints on the fraction of 
the Galactic Halo mass due to compact objects in the 
mass range [10 -2 , 1] Mq. These limits, along with the fact 
that the two SMC events observed so far are probably due 
to SMC deflectors, suggest that lmc and SMC self-lensing 
may dominate the event rate. 

Key words: : Galaxy: halo - Galaxy: kinematics and dy- 
namics - Galaxy: stellar content - Magellanic Clouds - 
dark matter - gravitational lensing 



1. Introduction 



In 1986, B. Paczyhski ( Paczyhski 1986 ) proposed to use 
gravitational microlensing as a tool to probe the dark mat- 
ter of the Galactic Halo. A few years later, several experi- 
ments started looking for microlensing events towards the 
Large Magellanic Cloud (lmc). Several candidates have 
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now b een observed by t h e MACHO and EROS collabora- 
tions ( Alcock et aL 1995 , Aubourg et al. 1993 , Alcock et 
al. 1997a, |Ansari et al. 1996|) . 

The first results towards the LMC suggested an opti- 
cal depth of order half that required to account for the 
dynamical mass of the dark halo, and the typical Ein- 
stein radius crossing time of the events implied a mass of 
about 0.5 Mq for the lenses ( Alcock et al. 1997a ). The re- 
sults of the EROS photographic plates microlensing survey 
pointed to a somewhat lower optical depth (Ansari et al. 
1996). Such objects cannot be main-sequence stars, whose 
observed density is two orders of magnitude too low to 
account for the measured optical depth. Many hypotheses 
have been debated regarding the nature and location of 
these lenses, as LMC se lf-lensing QSahu 1994 |Wu 1994| ), 
tidal debris (Zhao 1998), intervening populations (Zarit- 
sky & Lin 1997), a warped Galactic Disk (Evans et al. 
1998), white dwarfs (Adams & Laughlin 1996, Graff et 
al. 1998), primordial black holes (Jedamzik 1998), all of 
which have brought forth numerous counter-arguments 
(Gould 1995, Alcock et al. 1997b, Gould 1998, Bennett 
1998, Beaulieu & Sackett 1998). No fully satisfactory pos- 
sibility has emerged, the main two trends being 0.5 Mq 
dark objects in the halo, such as primordial black holes, 
or 0.1 M Q main sequence stars near the LMC itself. The 
cosmological implications of these two hypotheses differ 
considerably. 

The Small Magellanic Cloud (smc) is highly valuable 
as a second target for halo microlensing searches. If the 
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lenses belong to the halo of our Galaxy, the optical depth 
and the distribution of event durations should be simi- 
lar to the ones observed towards the lmc. In the case of 
SMC or lmc self-lensing, however, the different dynamical 
properties of the two Clouds will give rise to differences in 
optical depth and event duration. 

Eros, whose setup was upgraded in 1996, is engaged in 
observations towards the lmc and the SMC. The analysis 
of the first year of observations towards the SMC yielded 
one event (Palanque-Delabrouille et al. 1998, Alcock et al 
1997c). A second event, whose source star was too faint to 
be in the EROS reference catalog, was alerted on by MACHO 
and actively monitored by all microlensing collaborations. 
It was generated by a binary deflector most probably lo- 



4 points above 4>q + la (below (f>o — la). The significance 
LP of a given variation is defined as: 



cated in the SMC itself (Afonso et al. 1998, Albrow et al. 
1999, Alcock et al. 1999, Rhie et al. 1999, Udalski et al. 
1999). This light curve, however, is not counted as an event 
in the following analysis which only uses stars identified 
on our template images, for which the efficiency can be 
estimated. 

We present here the details of the EROS 2 two-year 
analysis of data towards the SMC and discuss its implica- 
tions for the nature of the halo dark matter. 



2. Experimental setup and observations 

The telescope, camera, telescope operations and data re- 
duction are as described in Palanque-Delabrouille et al. 
1998 and references therein. Ten square degrees are mon- 
itored on the SMC. Note that since one CCD of the red 
camera was non functional, all the analysis is done on 
only 7 CCD's (i.e. a total field of 8.6 square degrees). 

The two-year data set contains 5.3 million light curves 
covering the period from July 1996 to March 1998. 



3. Data analysis 

The analysis of the two-year data set is similar to that 
of the first year (Palanque-Delabrouille et al. 1998). The 
major difference comes from the addition of rejection cuts 
requiring a minimum time coverage of the event (not pos- 
sible with the previous one- year time span) , and from the 
tuning of some cuts applied to parameters with a clear 
field-dependent distribution. The criteria are sufficiently 
loose not to reject light curves deformed by blending or 
by the finite size of the source, or events involving multi- 
ple lenses or sources. Most variable stars are rejected by 
at least two distinct cuts. 

As in the first year analysis, we define a positive (neg- 
ative) fluctuation as a series of data points that (i) starts 
by one point deviating by at least ler from the baseline 
flux (j)Q , (ii) stops with at least three consecutive points 
below <po + lcr (above 0o — If) and (iii) contains at least 



i=N 

LP = -J2 log 



Erfc 



V2j 



(1) 



where Xi is the deviation in cr's of the point taken at time 
ti and N the number of points within the fluctuation. We 
order the fluctuations along a light curve by decreasing 
significance. The cuts are described hereafter: 

Selection of microlensing candidates: 

— la: The main fluctuation detected in the red and blue 
light curves should be positive and occur simultane- 
ously in the two colors: if I is the time interval during 
which the data are more than la away from the base- 
line flux, we require (I rod n Iblue)/ (hed U Ibiue) > 0.2. 

— lb: To reject flat light curves with only statistical fluc- 
tuations, we require that on a given light curve 
LP(2 nd most significant fluct.) / LP(main fluct.) < 0.6 
in both colors .0 

— lc: We require that LP(main fluct.) > 30 in both 
colors. 

Rejection of variable stars: 

— 2a: To exclude short period variable stars that exhibit 
irregular light curves, we require that the RMS of the 
distribution of the deviation, in cr's, of each flux mea- 
surement from the linear interpolation between its two 
neighboring data points be small, typically less than 
2.2 (the exact value varying from field to field, with a 
looser limit set on outer fields). 

— 2b: To exclude variable stars which exhibit correlated 
fluctuations between the two colors, we calculate the 
correlation coefficient p between the "red" and the 
"blue" light curves, excluding points in the main fluc- 
tuation (enlarged by a 25% time margin on either side) 
so as to consider only the un-amplified part of the light 
curve. The cut is applied to 



cor = y/Np-2 x p/y/l-p? 

whose distribution is less sensitive than p to the num- 
ber of points N p . We require cor < 6. 

— 2c: We remove two under-populated regions of the 
color-magnitude diagram that contain a large fraction 
of variable stars (upper main sequence and bright red 
giants, see details in first year analysis). 

S/N improvement of the set of selected candidates: 

— 3a: We remove events with low signal-to-noise ratio by 
requiring a significant improvement of a microlensing 
fit (ml) over a constant flux fit (est), i.e. that 



1 This criterion has been relaxed slightly compared to the 
first year analysis. 
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2 

min 



[ X 2 (cst) - X 2 (ml)]/[x 2 (ml)/d.o.f.] > A X 
where d.o.f. is the number of degrees of freedom. As 
for cut 2a, the exact value of Ax m j n varies from field 
to field because of the variations in time sampling. 
Typically, Ax min ~ 350, with a higher limit on well- 
sampled fields. 

— 3b: We require that the maximum magnification in 
the microlensing fit be greater than 1.40. 

Time coverage of the event: 

— 4a: We require that the fitted time of maximum mag- 
nification be contained in the period of observation. 

— 4b: We require that the fitted value of the Einstein 
radius crossing time At < 300 days. 

Physical blending: 

— 5: If blending significantly improves the microlensing 
fit, then the fitted blending should be physical, i.e.: if 
X 2 (nobl.) — x 2 (bl-) > 10, we require that the fitted 
blending fraction /bi//tot > —0.3, in both colors. 

The tuning of each cut and the estimation of the effi- 
ciency of the analysis (with the correction due to blend- 
ing) is done with Monte Carlo simulated light curves, as 
described in Palanque-Delabrouille et al. 1998. The mi- 
crolensing parameters are drawn uniformly in the fol- 
lowing intervals: time of maximum magnification to G 
[^first — 350, ti as t + 350] days, impact parameter normal- 
ized to the Einstein radius uq G [0, 2] and time-scale 
At G [0, 350] days. Cuts 2a and 3a are tuned individually 
on each field so as to reject respectively ~ 15% and ~ 25% 
of the remaining Monte Carlo light curves. The impact of 
each cut on data and simulated events is summarized in 
table 1. 



Cut description 


Number of 

stars 
remaining 


Fraction of remaining 
stars removed by cut 


Data 


Simulation 


Stars analyzed 


5,307,774 






la: Simultaneity 


93,900 


98% 


76% 


lb: Uniqueness 


29,742 


68% 


11% 


lc: Significance 


9,978 


66% 


06% 


2a: Stability 


4,365 


56% 


16% 


2b: Correlation 


3,502 


20% 


02% 


2c: HR diagram 


2,996 


14% 


04% 


3a: Microlensing fit 


102 


97% 


24% 


3b: Magnification 


37 


63% 


16% 


4a: Time inclusion 


28 


24% 


20% 


4b: Duration 


22 


21% 


14% 


5: Blending 


1 


95% 


12% 



Table 1. Impact of each cut on data and simulated mi- 
crolensing events (allowing for blending) . Each fraction for 
cut n refers to the light curves remaining after cut (n — 1). 

There is one remaining candidate, the one already dis- 
cussed in the analysis of the first year data set. Its updated 
light curve is shown in figure [l]. 
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Fig. 1. Light curve of microlensing event EROS-SMC- 
97/1 (the fit shown does not include the periodic mod- 
ulation) . 



The efficiency of the analysis for events with an impact 
parameter uo < 1 and normalized to an observing period 
T Q bs of one year is summarized in table g, as a function of 
the Einstein radius crossing times At (in days). The main 
source of systematics could come from a biased estimation 
of the blending effect. We have studied various blending 
models and found the relative error to be less than 10%. 



At 


7 


22 


37 


52 


67 


100 


150 


250 


300 


350 


e 


4 


12 


16 


18 


21 


24 


25 


23 


13 


3 



Table 2. Efficiency (in %) of the analysis as a function 
of the time-scale At (in days) for events with uo < 1, 
normalized to T Q b s = 1 yr. We monitor N ^ s — 5.3 x 10 6 
stars. 



4. Limits on the contribution of dark compact 
objects to the Halo 

We fit a microlensing curve to the data of the candidate 
(see table ||), allowing for a periodic modulation of the 
source star as evidenced in the first year analysis. Blend- 
ing and time-scale being degenerate quantities whose 
common fit is quite sensitive to systematics, we set the 
blending to that measured by OGLE ( Udalski et al. 1997 ): 
/bl//tot ~ 24%, considering it to be a lower limit on the ac- 
tual contribution (due to possible additional non-resolved 
companions). This value of the blending fraction is com- 
patible with the absence of distorsion in the spectrum 



(3ahu 1998) 



uo 


to 


At 


Period 


^mod R 


^mod B 


0.424 
±0.004 


2568.8 
±0.8 


129 

±2 


5.126 
±0.002 


0.031 
±0.003 


0.022 
±0.003 



Table 3. Result of microlensing fit to the SMC event, with 
blending set to fu/f tot = 25%. to is the time of maximum 
magnification in days since Jan, 1 1990, At; the Einstein 
radius crossing time, in days. The \ 2 is 261 for 279 d.o.f. 
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In order to set limits on the contribution of dark ob- 
jects to the Halo, a Halo model must be used to obtain, for 
a given deflector mass, both a number of expected events 
and a distribution of event durations. We use the so-called 
"standard" halo model described in Palanque-Delabrouille 
et al. 1998 as model 1, and take into account the efficiency 
of the analysis given in the previous section. 

Moreover, to set conservative limits, we have assumed 
that the observed event is a halo event, in spite of the 
fact that it is most probably an SMC self-lensing event 
(Palanque-Delabrouille et al. 1998). 

Assuming a standard halo model with a mass fraction 
/ composed of dark compact objects having a single mass 
M, the likelihood of observing at most one event, with a 
duration less probable than the observed one (Ai b s ), is 



L(M) = e- f " M (l + /JVmPj* (Afobe)) 



(2) 



with 



A/, 



V M (At ohs )= P M (At)dAt + P M {At)dAt (3) 

JO JAt 2 

Aii and Ai2 are the two durations with the same proba- 
bility as Ai bs, defined by: 



P M (Aii) = P M (Ai 2 ), Aii < Ai 2 , 
(Aii = Ai obs or Ai 2 = Ai obs ) 



(4) 



Pw(Ai) is the distribution of expected event durations, 
taking efficiencies into account, and Nm is the expected 
number of events.^] 

Figure || shows the 95 % exclusion limit derived from 
this likelihood function on /, the halo mass fraction, at any 
given mass M, i.e. assuming all deflectors in the halo have 
mass M. For comparison, the exclusion curve obtained 
considering that the sole event detected does not belong 
to the halo (thin dashed line on the plot) is also shown. 
The limit is very similar to the one event one, since the 
very long duration of the event pushes it towards large 
masses, where our efficiency starts to drop. 
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Fig. 2. Exclusion diagram at 95 % CL. for the standard 
halo model (4 x 10 11 M Q inside 50 kpc). The dashed line 
is the limit from EROS 1 lmc and SMC data (Renault et 
al., 1997), the solid line is the limit from the SMC data 
described in this paper. The MACHO 95 % CL. accepted 
region is the shaded area, with the preferred value indi- 
cated by the cross (Alcock et al. 1997a). The thin dashed 
line corresponds to the limit obtained assuming we ob- 
served no halo events. 

be obtained on the deflector location — through paral- 
lax for this 97-SMC-l event (Palanque-Delabrouille et al. 
1998), and deflector binarity for 98-SMC-l and LMC-9 
(but see the discussion in Bennett et al. 1996) — seem 
to be produced by deflectors located in the Clouds them- 
selves. This suggests that self-lensing may be the domi- 
nant source of the observed events. We are continuing to 
accumulate data towards the LMC and SMC in order to 
determine definitively the location of the deflectors. 
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5. Discussion and conclusion 

The analysis of two years 



lysis ot two years of EROS 2 SMC data has yielded 
a single microlensing event. This allows us to put new 
nts on the traction ot the halo made ot obiccts m 



constraints on the traction ot the halo made ot objects m 
the range [iu _2 Mq, IM0J, excluding m particufar at the 
95 % CL. that more than 50 % of the standard halo be 
made of 0.5 M Q objects. 

One should also note that all the microlensing events 
towards the LMC and the SMC for which information could 
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